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Primary gas thermometry by means of laser absorption spectroscopy: 
Determination of the Boltzmann constant 
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We report on a new optical implementation of primary gas thermometry based on laser absorption 
spectrometry in the near infrared. The method consists in retrieving the Doppler broadening from 
highly accurate observations of the line shape of the R(12) v\ -\- 2 V2^ + ^'a transition in CO2 
gas at thermodynamic equilibrium. Doppler width measurements as a function of gas temperature, 
ranging between the triple point of water and the gallium melting point, allowed for a spectroscopic 
determination of the Boltzmann constant with a relative accuracy of ^LGxlO"*. 

PACS numbers: 06.20.Jr, 33.20.-t, 42.62.Fi 



Atomic and molecular spectroscopy have played a ma- 
jor role in fundamental metrology since the advent of 
highly monochromatic, tunable laser sources. Spectro- 
scopic experiments with a metrological goal mostly con- 
sist in precision measurements of absolute frequencies, 
such as the center frequency of the hydrogen ls-2s line 
[l|, I4I and the strontium intercombination line [3|, just to 
mention a few significant examples. 

Recently, high resolution laser absorption spectroscopy 
has been proposed for a direct determination of the Boltz- 
mann constant [J]. In this case, a frequency interval has 
to be measured very accurately. In fact, the method re- 
lies on a determination of the line width of a Doppler- 
broadened absorption profile for an atomic or molecular 
line in a rarefied gas at thermodynamic equilibrium. The 
current value of fee, 1.3806504 (24) x lO^^a jj^-i^ ^^^_ 
commended by the Committee on Data for Science and 
Technology (CODATA) [3| , comes from the ratio between 
the molar gas constant R and the Avogadro number A^a 
[y]. Its relative uncertainty, 1.7 x 10~^, is mostly due 
to the uncertainty on R, whose accepted value is essen- 
tially given by that obtained in 1988 by Moldover et al. 
[7| through a fascinating experiment based on acoustic 
gas thermometry in argon. There is presently a strong 
interest in new and more direct determinations of fee in 
view of a possible new definition of the unit kelvin [8| . In 
the current International System of Units (SI), the unit 
kelvin is linked to an artefact, which is far from being 
invariant in space and time. As already done for the 
unit metre by fixing the speed of light c, the unit kelvin 
could be redefined through the thermal energy by fixing 
the Boltzmann constant. This new approach would make 
the definition more general and independent of any ma- 
terial substance, technique of implementation, and tem- 
perature. Nonetheless, before fixing its value, fee must be 
determined with an improved uncertainty, possibly using 
different methods. Presently, the most accurate methods 
to measure k^ are based on acoustic gas thermometry, ex- 
ploiting the relation between the speed of sound in a gas 
and the thermodynamic temperature, and on dielectric 
constant gas thermometry, based on the virial expansion 



of the Clausius-Mossotti equation and applicable only to 
helium, for which QED provides the static electric dipole 
polarizability [9|. Nonetheless, the experimental inves- 
tigation of new primary thermometric methods, based 
on radically different principles, is strongly encouraged 
by the Consultative Committee for Thermometry of the 
International Committee for Weights and Measures, as 
reported in its Recommendation T 2 of 2005. 

Doppler-broadened laser absorption spectroscopy can 
be implemented in such a way as to act as a primary 
thermometric method, with the advantage of being con- 
ceptually simple, applicable to any gas at any tempera- 
ture, in whatever spectral region. In contrast to other 
methods also based on electromagnetic radiation mea- 
surements, such as total radiation thermometry [10| . this 
new approach does not require absolute radiation deter- 
minations. It is well known that the Doppler width, Avjj 
(FWHM), of a line (with a center frequency i^) in an 
absorbing molecular gas at thermodynamic equilibrium 
depends on the temperature T through the equation: 
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where m is the mass of the molecule. 

The first spectroscopic determination of fee has been 
recently performed in the mid-infrared by the Daussy et 
al. on the 1^2 asQ(6,3) rovibrational line of the ammo- 
nia molecule "NH3 at a frequency of 28 953 694 MHz 
llj . The absorption profile was observed with a maxi- 
mum signal-to- noise ratio (S/N) of 10^, in the pressure 
range between 1 and 10 Pa, at a temperature of 273.15 
K, using a CO2 laser frequency stabilized on a OSO4 line. 
Under these low-pressure conditions, which are close to 
the Doppler limit, by measuring the width of the absorp- 
tion line as a function of the pressure and extrapolating 
to zero pressure, it was possible to deduce the Doppler 
width, which yielded an evaluation of the Boltzmann con- 
stant with a relative uncertainty of 1.9x10""' [ll|. This 
approach allows for a very simple spectral analysis but 
requires an accurate determination of the absolute pres- 
sures. 




Vacuum pump PG 

Figure 1: A sketch of the intensity-stabihzed laser absorption 
spectrometer. ECDL: extended-cavity diode laser; M: mirror; 
BS: beam splitter; FP: fiber port; Ph: photodetector; PG: 
pressure gauge. The cell is closed by an oil-free solenoid valve, 
placed inside the chamber. 



Here we report on a radically different implementation 
of laser absorption spectroscopy for primary gas ther- 
mometry. In particular, we demonstrate that it is possi- 
ble to retrieve the gas temperature from a molecular ab- 
sorption profile even when the gas pressure is sufficiently 
high that the line shape is far from the Doppler limit, 
but sufficiently small that one can neglect the averag- 
ing effect of velocity-changing collisions, so that the line 
shape is given by the exponential of a Voigt convolution. 
We perform absorption spectroscopy in the near-infrared 
on a CO2 gas sample at thermodynamic equilibrium us- 
ing a distributed feed-back (DFB) diode laser, probing 
the R(12) component of the vi + 2V2' + v-i combination 
band. In contrast to the NH3 molecule, our molecular 
target does not exhibit any hyperfine structure. 

The experimental setup is depicted in Fig. [TJ The 
diode laser, emitting on a single mode at 2.006 \jxa. with 
a maximum power of 3 mW, is mounted in a mirror 
extended-cavity configuration. In this scheme, a par- 
tially reflecting mirror, with a 50% reflection coating on 
one side and an antirefiection (AR) coating on the other, 
provides the optical feedback for line narrowing down 
to '^1 MHz. As will be shown later on, this is suffi- 
cient for the aim of the present experiment. CoUimated 
by a double-aspheric lens, the laser beam first passes 
through an acousto-optic modulator (AOM) and then 
is injected into a single-mode, polarization-maintaining 
optical fiber. The latter is employed as a high-quality 
spatial filter to produce an optimum TEMqo transverse 
mode. The AOM is used as an actuator within a servo 
loop in order to keep the intensity of the laser beam con- 
stant over a laser frequency scan as wide as 1 wavcnum- 
ber. Driven by a 1 W radiofrequency signal at 80 MHz, 
the AOM deflects about 70% of the laser power from 
the primary beam to the first diffracted order, which is 
used for the experiment. The intensity control feedback 
loop ensures a stability level better than 10^"* (measured 
value over a laser frequency span of about 0.5 cm~^). 
The servo loop is active within an electronic bandwidth 



of 50 kHz. After being collimated again, the laser beam 
is divided into four parts by a series of beam splitters. 
The first beam is sent to a high-finesse, confocal optical 
resonator that provides accurate frequency markers for 
calibration purposes. The second goes to a 5 cm tick ger- 
manium etalon, which allows for a coarse check of the fre- 
quency calibration. The third beam impinges on a moni- 
tor photodiode, whose output signal is used as the input 
signal to our servo to control the rf power driving the 
AOM. Finally, the remaining portion of the laser beam 
passes through a homemade, temperature-stabilized sam- 
ple cell, closed at each end by an AR-coatcd BK7 win- 
dow. Temperature-stabilized, extended-wavelength, In- 
GaAs photodctectors are employed to monitor the laser 
radiation. The optical resonator consists of a pair of 
identical, high-quality, dielectric mirrors, with a radius 
of curvature of 50 cm and a reflectivity of 99.5%. The 
mirror mounts, tightly locked over a stainless steel base, 
fit into the two ends of a glass tube, which is periodically 
evacuated. The absorption cell is 10.5 cm long and is 
housed inside a stainless steel vacuum chamber, which 
is equipped with a pair of optical windows of the same 
type as those of the cell. Consisting of a cylindrical cav- 
ity inside an aluminium block, with inner and external 
surfaces carefully polished, the cell is temperature stabi- 
lized by means of four Peltier elements. Three precision 
platinum resistance thermometers (PtlOO) measure the 
temperature of the cell's body while a proportional inte- 
gral derivative controller is used to keep the temperature 
uniform along the cell and constant within 40 mK, over a 
time interval of ~2 hours. This active system also allows 
us to vary the gas temperature between 270 and 330 K. 

The PtlOO thermometers were calibrated at the triple 
point of water and at the gallium melting point with 
an overall accuracy better than 0.01 K. Both during the 
calibration and when placed in the absorption cell, the 
thermometers are fed by a 1 mA current. The correction 
due to the self heating effect has been applied together 
with the associated uncertainty. 

The sample cell is filled with CO2 gas (with a nom- 
inal purity of 99.999%) at a pressure between 70 and 
130 Pa, measured using a 1300 Pa full-scale capacitance 
gauge with a 0.25% accuracy. A turbomolccular pump 
is used to periodically evacuate the sample cell and cre- 
ate high-purity conditions. Data acquisition is performed 
by means of a digital oscilloscope (Tektronix TDS7104), 
with a 13-bit vertical resolution (on the voltage range 
covered by the absorption profile) and a total number of 
points for each spectrum equal to 5000. Spectral averag- 
ing over 25 consecutive scans is performed, the scan rate 
being 5 Hz. The resulting spectrum is then transferred 
to a personal computer. Simultaneously, the transmis- 
sions from both the optical resonator and the etalon are 
recorded. 

An example spectrum is shown in Fig. [51 The observed 
S/N was 1600:1, with an equivalent noise detection band- 
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Figure 2: Transmission spectrum for the R(12) line at ~100 
Pa and 275.12 K. The absolute residuals are also shown, as a 
result of a non-linear least-squares fit to a Voigt convolution. 
The retrieved Doppler width was 268.01 MHz. 
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Figure 3: Repeated spectroscopic measurements of the abso- 
lute temperature and a histogram of the determinations. The 
mean value is 298.77 K, the standard deviation 0.89 K, and 
the standard error 0.13 K. 



width of 600 Hz and a power incident upon the cell of 
50 /.tW. A MATLAB code is used to carry out the line 
shape fitting and retrieve the Doppler width in a fully au- 
tomated way for any set of repeated acquisitions. Partic- 
ular attention is paid to the frequency calibration of the 
spectra. In fact, prior to each fit, the same code analyzes 
the transmission spectrum from the confocal resonator 
in order to construct a linear frequency scale that has to 
be matched to the corresponding absorption spectrum. 
This is done by means of a cubic spline interpolation of 
the peak frequencies versus their index in the data ar- 
ray. Obviously, an accurate determination of the cavity 
free-spectral-range, /fsR: is also necessary. For that pur- 
pose, two different methods were implemented, whose 
results were in total agreement. The first one was based 
on the simultaneous detection of absorption on both the 
zero- and first-order beams as provided by the AOM, 
once the former was sent to a second absorption cell (see 
Fig. [1]). The two cells were filled with CO2 gas at the 
same pressure and at room temperature. Because of the 
frequency shift between the two beams (80517. 42±0. 12 
kHz, as measured by a universal counter), the center fre- 
quencies of the two line profiles, simultaneously recorded 
along with the transmission comb from the optical cav- 
ity, are separated by an exactly known quantity. The 
second method consisted in making a broad laser fre- 
quency scan (~0.5 cm~^) in order to observe a pair of 
CO2 absorption lines, namely the R(8) 1^1 4- 21^2" + '^3 
and R(27) 1^1 + 2>Vr^ — v^ + vz transitions, whose cen- 
ter frequencies are accurately known. These frequencies, 
as well as that of the R(12) line, were taken from Ref. 
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Thus, as a result of a weighted mean between the 
two values, /fsr was found to be 149.56±0.01 MHz. In 
our operation conditions, the effects of temperature vari- 
ations and changes in the cavity alignment on /fsr could 
be neglected. As it will be shown later on, the relative 



uncertainty on the cavity free-spectral-range, multiplied 
by a factor of 2, will contribute to the overall accuracy 
of the spectroscopic determination of /cb . The linearity 
of the frequency scale was periodically checked using the 
transmission spectra from the germanium etalon. Devi- 
ations from linearity, which could affect the symmetry 
of the absorption profiles, were found to be negligible 
at the present experimental accuracy. Coming back to 
the fit of the absorption line shape, the normalized Voigt 
function was computed using the approximation of Wei- 
deman [l3|]. It was confirmed that this approximation is 
equal to the Voigt convolution, as tabulated in Ref. [ij, 
to better than 1 part in 10^. The nonlinear least-squares 
fitting code was based on a Levenberg-Marquardt algo- 
rithm. The code was tested on simulated spectra, as- 
suming a gas pressure of 130 Pa at a temperature of 
296 K, and adding a random noise to the computed line 
shape. As one could expect, the relative accuracy in the 
retrieval of the Doppler width strongly depends on the 
signal-to-noise ratio. In fact, relative deviations between 
the real value and the retrieved one of 2.6xl0~^ and 
5.2xl0~^ (mean values over 10 simulations) were found 
for S/N=1600 and 1.6x10^, respectively. 

The exponential of a Voigt function is able to fit 
our measured spectra within the experimental noise, as 
shown in Fig. [51 Reported in the bottom part of the fig- 
ure, the residuals evidence that the noise is dominated 
by periodic fluctuations due to spurious etalon effects, 
whose root-mean-square value is of the order 10~^ V. 
Owing to our intensity stabilization procedure, the back- 
ground signal for the Voigt fit could be optimally repro- 
duced by simply adding a flat baseline to the line shape 
function. Apart from the Doppler line width, the re- 
maining free parameters were the line center frequency, 
the homogeneous width and an amplitude factor, the lat- 
ter representing the product of the line strength factor, 
the gas pressure and the absorption path length. In our 




Figure 4: A plot of the quantity aAu^ as a function of the gas 
temperature, where a = mc /{8u In 2) and aAi/j^ represents 
the thermal energy. The slope of the weighted best-fit line 
directly provides the value of ks- Each point represents the 
mean value over a number of repeated determinations, the 
error bars being the standard error (3cr) . These data resulted 
from more than 1500 spectra, acquired over five consecutive 
weeks. 



experimental configuration, saturation and transit time 
broadening could be neglected. It is worth noting that 
the scan goes quite far into the line wings, as demon- 
strated a posteriori by the area under the normalized 
fitting function, which was found to be larger than 0.99. 

In Fig. [21 an example of spectroscopic determinations 
of the gas temperature, over a time interval of ~4 min- 
utes, is reported. The histogram, on the right side of 
Fig. [3l essentially shows that the distribution is gaussian. 
Simultaneously to the spectral acquisitions, the PtlOO 
temperatures were also recorded, providing a mean value 
of 298.68 K, with a standard deviation of '^10 mK. With 
very few exceptions, a good agreement was found be- 
tween the spectroscopic and PtlOO temperatures, typical 
deviations being much smaller than 1 K. 

Doppler width determinations were repeated as a fun- 
ction of the gas temperature, in the range between the 
triple point of water and the gallium melting point. This 
approach was mainly adopted to cancel out possible sys- 
tematic errors arising from residual etalon effects due to 
the cell's windows. In Fig. [U the experimental values of 
the quantity gj^^Ai/^ and the corresponding temper- 
atures, as provided by the PtlOO thermometers, are plot- 
ted. The molecular mass was calculated from the relative 
atomic masses of ^^O and ^^C, the molar mass constant, 
and the Avogadro number [o] . A linear fit of these data 
directly provides the spectroscopic value for the Boltz- 
mann constant, namely /cb = 1.38058(22) x 10"^^ JK"^ 
The relative uncertainty in the temperature values, also 
including the short-term fluctuations, were found to be 
a factor of 10 smaller compared to the relative error of 
a typical spectroscopic determination of the temperature 
itself. For this reason, there errors were ignored in the 



linear fit. The 1-a uncertainty on fcs includes the statis- 
tical error on the slope of the best-fit line as well as the 
uncertainty on the frequency scale. Our determination 
agrees with the CODATA value. 

In conclusion, we have demonstrated that highly accu- 
rate, intensity-stabilized, laser absorption spectrometry 
in the near-infrared represents an extremely powerful tool 
for primary gas thermometry. Our experiment allowed 
us to provide a spectroscopic measurement of the Boltz- 
mann constant, with a relative accuracy of 1.6 xlO"'*. 
The absorption features of our molecular target around 
2-fj.m were found to be a very good choice mainly be- 
cause of the insensitivity of the Doppler width to the gas 
pressure. More precisely, for CO2 pressures below a few 
hundred Pa, we did not observe any influence from Dicke 
narrowing or speed-dependent effects. This is true at the 
present precision level. Once, in the near future, we will 
be able to retrieve the Doppler width with a much smaller 
uncertainty, it will be necessary to confirm its behavior 
as a function of the gas pressure. 

Future developments include the use of a pair of phase- 
locked extended-cavity diode lasers to improve our capa- 
bility of measuring laser frequency variations, as well as 
the implementation of a new detection method based on 
fast amplitude modulation of the laser beam and phase- 
sensitive detection to enhance the signal-to-noise ratio. 
We are also building a new temperature-stabilized cell, 
which should ensure a constant and homogeneous gas 
temperature within 1 mK. After these significant im- 
provements, our method could contribute to the new def- 
inition of the kelvin, in the near future. In fact, we will no 
longer be limited by the accuracy of the frequency scale. 
On the other hand, the S/N is expected to increase by 
a factor of 25, provided that spurious etalon effects are 
eliminated. We also plan to exploit the excellent linearity 
of standard InGaAs detectors by shifting the operation 
wavelength down to 1.4 /im. According to our simula- 
tions, it should be possible to reach an accuracy level of 
a few parts in 10^. 
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